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ABSTRACT 

We present the first analysis of extended stellar kinematics of elliptical galaxies where a Yukawa- 
like correction to the Newtonian gravitational potential derived from /(i?)-gravity is considered as 
an alternative to dark matter. In this framework, we model long-slit data and planetary nebulae data 
out to 7 R e g of three galaxies with either decreasing or flat dispersion profiles. We use the corrected 
Newtonian potential in a dispersion-kurtosis Jeans analysis to account for the mass-anisotropy degen- 
eracy. We find that these modified potentials are able to fit nicely all three elliptical galaxies and the 
anisotropy distribution is consistent with that estimated if a dark halo is considered. The parameter 
which measures the "strength" of the Yukawa-like correction is, on average, smaller than the one 
found previously in spiral galaxies and correlates both with the scale length of the Yukawa-like term 
and the orbital anisotropy. 

Subject headings: galaxies : kinematics and dynamics - galaxies : general - galaxies : elliptical and 
lenticular - cosmology: theory 



1. INTRODUCTION 

The "concordance" ACDM cosmological model, which 
includes some unseen Cold Dark Matter (DM) and a cos- 
mological constant (A) acting as a repulsive form of Dark 
Energy (DE) , has been remarkably successful in explain- 
ing the formation and evolution of cosmological struc- 
tures at different scales (e.g., Springel et al. 2006). 

However, at cosmological scales, the cosmological con- 
stant as "vacuum state" of the gravitational field is about 
120 orders of magnitude smaller than the value predicted 
by any quantum gravity theory (Weinberg 1989) and 
comparable to the matter density (coincidence problem), 
even if they evolved decoupled in the history of the uni- 
verse. 

In addition, looking at the galaxy scales there are a 
few critical issues yet to be solved, which are giving hard 
time to the whole ACDM framework. 

Since the discovery of the flat rotation curves of spiral 
systems, galaxies have been the most critical laboratory 
to investigate the gravitational effects of the DM halos, 
to be compared against the expectation of the cosmo- 
logical simulations (Navarro et al. 1997, NFW hereafter, 
Burkert 1995, Navarro et al. 2010). Here, the ACDM 
model is not able to fully explain the shallow central 
density profile of spiral and dwarf galaxies (Gilmore et 
al. 2007; Salucci et al. 2007; Kuzio de Naray et al. 2008). 
Early-type galaxies (ETGs hereafter) have been proven 
only recently to be consistent with ACDM predictions 
(and WMAP5 cosmological parameters, e.g. Komatsu 
et al. 2009), from their centers (Napolitano et al. 2010) 
to their peripheries (Napolitano et al. 2011, N+ll here- 
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after), although there are also diverging results showing 
that ETGs in some cases have too high (Buote et al. 
2007) or too low (e.g. Mandelbaum et al. 2008) concen- 
trations. 

This very uncertain context has been a fertile soil for 
alternative approaches to the so-called "missing mass". 
The basic approach is that the Newtonian Theory of 
Gravity, which has been tested only in the Solar System, 
might be inaccurate on larger (galaxies and galaxy clus- 
ters) scales. The most popular theory investigated so far, 
the Modified Newtonian dynamics (MOND) proposed by 
Milgrom (1983), is based on phenomenological modifica- 
tions of Newton dynamics in order to explain the flat 
rotation curves of spiral galaxies, and passed a number 
of observational tests (Ferreira & Starkman 2009), in- 
cluded ETG kinematics (Milgrom & Sanders 2003; Tiret 
et al. 2007; Kroupa et al. 2010; Cardone et al. 2010; 
Richtler et al. 2011). Only lately it has been derived in 
a cosmological context (Bekenstein 2004) . 

A new approach, motivated from cosmology and 
quantum field theories on a curved space-time, has 
been proposed to study the gravitational interaction: 
the Extended Theories of Gravity (Capozziello 2002; 
Capozziello & Faraoni 2011). In particular, the so called 
/(i?)— gravity seem to have passed different observational 
tests like spiral galaxies' rotation curves, X-ray emis- 
sion of galaxy clusters and cosmic acceleration (see e.g. 
Capozziello et al. 2007a, C+07 hereafter, Capozziello et 
al. 2009, C+09 hereafter, Capozziello et al. 2008). This 
approach is based on a straightforward generalization 
of Einstein theory where the gravitational action (the 
Hilbert-Einstein action) is assumed to be linear in the 
Ricci curvature scalar R. In the case of /(i?)-gravity, 
one assumes a generic function / of the Ricci scalar R 
(in particular analytic functions) and asks for a theory 
of gravity having suitable behaviors at small and large 
scale lengths. 

As shown in Capozziello et al. (2009), analytic f(R)- 
models give rise, in general, to Yukawa-like corrections to 
the Newtonian potentials in the weak field limit approx- 
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imation (see also Lubini et al. 2011). The correction in- 
troduces a new gravitational scale, besides the standard 
Schwarzschild one, depending on the dynamical structure 
of the self-gravitating system. 

Here we want to test these Yukawa-like gravitational 
potentials against a sample of elliptical galaxies. This 
approach has been proposed earlier, in a phenomcno- 
logical scheme for anti-gravity, to model flat rotation 
curves of spiral galaxies (Sanders 1984), and recently, 
in f(R) theories to model disk galaxies combined with 
NFW haloes (see Cardone & Capozziello 2011). The test 
we are proposing at galaxy scales is crucial: reproduc- 
ing kinematics and then dynamics of these very different 
classes of astrophysical systems in the realm of the same 
paradigm is needed to test these new gravitational the- 
ory as an alternative to DM which has not been definitely 
found out at fundamental level. 

The layout of the paper is the following. In §2, we 
sketch the main ingredients of /(i?)-gravity deriving, in 
the weak field limit, the Yukawa-like corrected gravita- 
tional potential. §3 is devoted to the high-order Jeans 
analysis suitable for ellipticals. The dispersion-kurtosis 
fitting and the data sample are presented in §4. Discus- 
sion and conclusions are in §5. 

2. POST- NEWTONIAN POTENTIALS FROM 
F(R)-GRAVITY 

We are interested in testing a class of modified poten- 
tials which naturally arise in post-Newtonian approxi- 
mation of f(R)— gravity for which no particular choice 
of the Lagrangian has been provided. 

The starting point is a general gravity action of the 
form 



/ 



d 4 x. 



f(R) + XC r , 



(1) 



where f(R) is an analytic function of Ricci scalar, g is the 

„ 16ttG . , 
determinant oi the metric g„ u , X = — — is the gravita- 

c 4 

tional coupling constant, and C m describes the standard 
fluid-matter Lagrangian. Such an action is the straight- 
forward generalization of the Hilbert-Einstein action ob- 
tained as soon as f(R) = R- 

In Capozziello et al. (2009, and reference therein) it 
has been shown that if one solves the field equations in 
the weak field limit under the general assumption of an 
analytic Taylor expandable f(R) functions of the form 

f(R) ~ / + hR + hR 2 + fsR 3 + ... (2) 
the following gravitational potential arises 



$ ( GM | ^iW e ' f N ) 
V fir 6 r J 



(3) 
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where L = — — , fi and fa are the expansion co- 
Ji 

efficients obtained by Taylor expansion. We note that 
the L parameter is related to the effective mass m = 
(-3/L 2 )- 1 / 2 = (2/ 2 /A) 1/2 and can be interpreted also 
as an effective length. 

From Eq. 3, the standard Newton potential is recov- 
ered only in the particular case f(R) = R- Furthermore, 



the parameters f\ and fi and the function b\ represent 
the deviations with respect the standard Newton poten- 
tial. On the Solar system scale, it has been shown that 
Yukawa-like deviation from the pure Newtonian poten- 
tial are not in contradiction with classical tests of Gen- 
eral Relativity (see e.g. Capozziello & Tsujikawa 2008), 
thanks to the so-called Chameleon mechanism (Khoury 
& Weltman 2004). In particular, fx and fi parameters 
are expected to allow the regular Newtonian potential, 
while at larger scales they can assume non-trivial values 
(e.g. fi 7^ 1, S\(t) 7^ 0, £ 7^ 1, see Capozziello et al. 
2007b, Capozziello et al. 2009). 
Eq. (3) can be recast as 



$(r) 



GM 
(1 + 6)r 



(4) 



where the first term is the Newtonian-like part of the 
potential associated to baryonic point-like mass M/(l + 
5) (no DM) and the second term is a modification of the 
gravity including a "scale length", L associated to the 
above coefficient of the Taylor expansion. If S = the 
Newtonian potential is recovered. Comparing Eqs. 3 and 
4, we obtain that 1 + S = /1, and S is related to Si(t) 
through 

6GM S 

61 = Tp~YT5 (5) 

where 6GM/L 2 and <5i can be assumed quasi-constant. 
From Eq. 5, it turns out that L cx y/—S/(l + S). Due to 
the arbitrarity of Si(t), the actual value of the S param- 
eter can assume any values, however, in order to have a 
Yukawa potential with a non imaginary exponent (i.e. L 
must be real) it is required that £ < or — 1 < 6 < 0. As 
comparison, Sanders (1984) adopted the same potential 
as in Eq. 4 under the assumption of anti-gravity gen- 
erated by massive particles (of mass m ) carrying the 
additional gravitational force. In this case a typical scale 
length would naturally arise (L = h/niQC being a Comp- 
ton length) and a — 1 < S < would provide a repulsive 
term to the Newtonian-like term producing flat rotation 
curves at r ^> L as observed in spiral galaxies. In partic- 
ular, for a small sample of spiral systems Sanders (1984) 
found -0.95 £ 5 £ -0.92. 

Here we want to test the modified potential as in Eq. 
4 in elliptical galaxies and check whether it is able to 
provide a reasonable match to their kinematics and how 
the model parameters compare with the results obtained 
from spiral systems. We construct equilibrium models 
based on the solution of the radial Jeans equation (see §3) 
to interpret the kinematics of planetary nebulae (PNe, 
see Napolitano et al. 2002, 2005; Romanowsky et al. 
2003; Coccato et al 2009) which are the only stellar- 
like tracers for galaxy dynamics available in ETGs out 
to <~ 5 — 10 effective radii (R e s)- 

We will use the inner long slit data and the extended 
PN kinematics for three galaxies which have published 
dynamical analyses within DM halo framework: NGC 
3379 (Douglas et al. 2007; De Lorenzi et al. 2009, DL+09 
hereafter), NGC 4494 (Napolitano et al. 2009, N+09), 
NGC 4374 (N+ll). The decreasing velocity dispersion 
profiles of the first two galaxies have been modeled with 
an intermediate mass halo, logMvi r ~ 12 — 12.2M , 
with concentration c v i r = 6 — 8 and a fair amount of 
radial anisotropy in the outer regions. For NGC 4374, 
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Fig. 1. — Circular velocity produced by the modified potential 
in Eq. 4 for the two galaxies N4494 (top) and N4374 (bottom). 
In both cases the M/L„ has been fixed to some fiducial value (as 
expected from stellar population models and Kroupa 2001 IMF): 
Af/L. = 4.3T QiS for NGC 4494 and M/L, = 5.5T Qi v for NGC 
4374. The potential parameters adopted are: L = 250" and <5=0, 
-0.65, -0.8, -0.9 (lighter to darker solid lines) and L = 180" and S=- 
0.8 (dashed lines). The dotted line is a case with positive coefficient 
of the Yukawa-like term and L = 5000" which illustrates that 
positive <5 cannot produce flat circular velocity curves. Finally 
some reference Navarro-Frenk- White (NFW) models are shown as 
dot-dashed lines. 

having a rather flat dispersion profile, a more massive 
(adiabatically contracted) halo with logM v ; r <~ 13.4M 
and c v ir ~ 7 was required with a negligible amount of 
anisotropy in the outer regions. These models turned out 
to be in fair agreement with the expectation of WMAP5 
c v ; r — relation (N+ll) and with a Kroupa (2001) 
IMF, making this sample particularly suitable for a com- 
parison with an alternative theory of gravity with no-DM 
as we want to propose here. 

Before we go on with detailed stellar dynamics, we 
show in Fig. 1 the circular velocity of the modified po- 
tential as a function of the potential parameters L and 
5 for NGC 4494 and NGC 4374. As for the spiral galax- 
ies, negative values of the <5 parameter make the circular 
velocity more and more flat also reproducing the typical 
dip (e.g. NGC 4374) of the circular velocity found for 
the DM models (dot-dashed curves) of the most mas- 
sive systems. On the contrary, positive 5 values cannot 
produce flat circular velocity curves (see Fig. 1). 

2.1. A consistency check with galaxy scaling relations 

To conclude the inspection of the modified potential as 
in Eq. 4 here we want to show that, beside flat rotation 
curves, this also naturally accounts for fundamental scal- 
ing relations of galaxies: the Tully-Fishcr (TF) relation 
for spirals and Faber-Jackson (FJ) relation for ETGs. 

Both relations connect the total mass M of galaxies 
with some characteristic velocity defining the kinetic en- 



ergy of the systems (i.e. the maximum rotation velocity, 
for spirals and the central velocity dispersion, a , 
for ETGs). In cither cases the kinematical quantities 
involved are proportional to the circular velocity of the 
systems through some "structure" constant, thus the ar- 
guments below apply to galaxies in general. 

Although the point-like version of the potential implies 
that the circular velocity v c scales with mass as M ~ 
w c,max ( as pointed out by Sanders 1984), if one derives 
the circular velocity for an extended galaxy this can be 
generalized as 



(GM tot /r,) x f(r/r*;6,L/r*) 



(6) 



where r* is a characteristic radius (e.g. the disk length for 
spirals or the effective radius encircling half of the galaxy 
light for ETGs), /(r/r*; S, L/r*) is a generic function 
which includes the radial dependence of the enclosed 
mass and the above Yukawa-like term. This function 
is defined such as, for 6 = 0, it gives v^(r) = GM(r)/r 
as the usual Newtonian expression. It is easy to show 
that if galaxies are homologous the maximum of v c is 
reached at the same r/r*, for a given 5 and L/r* and 
this maximum can be written as 



<max = KM t , 



t/r* 



(7) 



where the constant K depends on the set of parameters 
{5, L, r*} adopted. In Eq. 7, though, M tot and 
are linked by the size-mass relation which is generally 
written as r* cx M t " t , from which Eq. 7 can be written 

as 

(8) 



Vmax <* M tot' 



The size-mass relation of spiral galaxies can be found in 
Persic et al. (1996, see also Thomas et al. 2009) to be 
r* cx M t ° ot 4 , while it is r* cx M t ° ot 6 for ETGs (e.g. Shen 
et al. 2003, Napolitano et al. 2005). This would give a 
TF slope of 3.33 and FJ slope of 5 which are both in the 
range of the observed relations (see e.g. McGaugh 2005 
and Nigoche-Netro et al. 2010 respectively) with the re- 
maining discrepancy being mainly due to the conversion 
factor to the observed quantities and non homologies. 

We finally remark that the TF relation has been found 
not to be conflicting with f(R) potentials in Capozziello 
et al. (2006), although the potentials from f(R) cx R n 
adopted there are just a series expansion of the Yukawa- 
like potential coming out from a more general polynomial 
f(R) as in Eq. 3. 

3. HIGH-ORDER JEANS ANALYSIS 

From the model point of view, the problem of fitting 
a modified potential as in Eq. 4 (which is formally self- 
consistent since the source of the potential is the only 
mass of the dynamical tracers, i.e. stars) implies the 
same kind of degeneracies between the anisotropy pa- 
rameter, j3 = 1 — <Tg/<Jr (where o~g and ay are the az- 
imuthal and radial dispersion components in spherical 
coordinates), and the non-Newtonian part of the poten- 
tial (characterized by two parameters like typical dark 
haloes) in a similar way of the classical mass-anisotropy 
degeneracy. We have shown (N+09, N+ll) that these 
degeneracies can be alleviated via higher-order Jeans 
equations including in the dynamical models both the 
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dispersion 6 (er p ) and the kurtosis (k) profiles of the trac- 
ers. 

In the following, we will use the assumption of spher- 
ical symmetry since galaxies in the sample are all E0- 
El for which, if one exclude the singular chance that 
they are all flattened systems seen face-on (see discus- 
sion in Sect. 8.1. of Douglas et al. 2007), the spher- 
ical approximation is good at 10% (Kronawitter et al. 
2000) 7 . Under spherical assumption, no-rotation and 
(3 = const (corresponding to the family of distribution 
functions f(E,L) = f a L~ 2/3 , see Lokas 2002 and refer- 
ences therein 8 , the 2-nd and 4-th moment radial equa- 
tions can be compactly written as: 



s(r) 



= r -2P 



x 2ft H(x)d2 



(9) 



100 200 300 100 200 300 100 200 

Radius [arcsec] Radius [arcsec] Radius [arcsec] 

Fig. 2. — Dispersion in km s _1 (top) and kurtosis (bottom) fit of the galaxy sample for the different f(R) parameter sets: the anisotropic 
solution (solid lines) is compared with the isotropic case (dashed line - for NGC 4374 and NGC 4494 this is almost indistinguishable from 
the anisotropic case). From the left, NGC 4494, NGC 3379 and NGC 4374 are shown with DM models as gray lines from N+09, DL+09 
(no kurtosis is provided), and N+ll respectively. 

This M/L+ — const might be a strong assumption 
to check further in a separate paper as it neglects the 
presence of stellar population gradients (see e.g. Tor- 
tora et al. 2010). However, colour (and M/L) gradients 
are generally stronger within i? e ff (see e.g. Tortora et al 
2010) and might mainly drive the best fit in the central 
regions, while they are possibly shallower outside (e.g., 
Tamura & Ohta 2003) where the f(R) parameters should 
be better constrained. In the following, j*(r) is derived 
by photometry presented in previous dynamical studies 
(i.e. DL+09, N+09, N+ll for NGC 3379, NGC 4494 and 
NGC 4374 respectively). 

Eqs. 9 are the ones interested by the potential mod- 
ification and include four free parameters to be best- 
fitted: the f(R) parameters {5, L}, the "dynamically in- 
ferred" stellar mass-to- light ratio Mj L* and the constant 
anisotropy f3 (see also §4). The solutions of Eqs. 9 on a 
regular grid in the parameter space are then projected to 
match the observed line-of-sight kinematical profile via 
ordinary Abel integrals (see N+09 for details). 

As mentioned earlier, Eqs. 9 are written under the as- 
sumption of a constant /3 with radius, which provides a 
average global anisotropy distribution over all the galaxy. 
As seen in previous analyses (e.g. N+09, DL+09 and 
N+ll), it is likely that this might not be a fair assump- 
tion, as [3 turns out to be constant somewhere in the 
outer regions, but strongly varying in the central radii. 
In this preliminary test we will skip this implementation 
of the models since we expect this to possibly improve 
the fit to the data in the central part only, where we do 
not expect the overall dynamics of being strongly ruled 
by the f(R) potential, whose parameters are the main 
focus of this work. Furthermore, we have shown previ- 
ously (see e.g. N+09 and N+ll) that the assumption of 
the constant or radial varying anisotropy did not strongly 
affect the determination of the other important parame- 
ter, the dynamically based stellar M/L. In the following 
we will take the f3 =const as fair estimate of the average 
galaxy anisotropy. 



where s(r) = {pa 2 ;pv*}, (3 is the anisotropy parame- 
ter, and 



dr dr 



respectively for the dispersion and kurtosis equations, 
being the latter n{r) — vf/crf. In the same equa- 
tions, is the spherical extended source version of 
the point-like potential as in Eq. 4 9 and p(r) is the 3D 
density of the tracer obtained by multiplying the depro- 
jection of the stellar surface brightness profile, j*(r), by 
some constant stellar mass-to-light ratio, M/L*. 

For the slow-rotating models we use the velocity, v TInB = 



\/v 2 + <t 2 as a measure of the velocity dispersion. 

7 The effect of non-spherical models is outside the scope of this 
paper, but details for NGC 3379 and NGC 4494 can be found in 
DL+09 and N+09. 

8 Here, there is the caveat that the solution of Jeans Equations 
does not ensure that the final distribution function is non negative 
and thus fully physical (see e.g., An & Evans 2006). 

9 This is obtained assuming the onion shell approximation: 
*(r) = /J" J 2 * <j>{r) r 2 sine dedipdr, see also Eq. 18 of C+09. 
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4. DISPERSION-KURTOSIS FITTING 

In Fig. 2 we show the dispersion and kurtosis profiles 
of the three galaxies with the f(R) models superimposed 
(solid lines). The fitting procedure is based on the simul- 
taneous \ 2 minimization of the dispersion and kurtosis 
profiles over a regular grid in the parameter space. The 
best-fit parameters are summarised in Table 1 together 
with some info of the galaxy sample. 

Overall the agreement of the model curves with data 
is remarkably good and it is comparable with models 
obtained with DM modeling (gray lines in Fig. 2). 

In all cases, the f(R) models allow to accommodate a 
constant orbital anisotropy /3 which is very close to the 
estimates from the DM models (see e.g. Table l 10 ). This 
is mainly guaranteed by the fit to the n{r) which does not 
respond much differently to the modified potential with 
respect the DM models. Thus, an important result of the 
analysis is that the orbital anisotropy is fairly stable to 
the change of the galaxy potential. In particular, the use 
of the kurtosis profiles has allowed us to solve the degen- 
eracy of the models and favor the anisotropic solutions 
for NGC 3379 and NGC 4494 (NGC 4374 being almost 
isotropic everywhere). Although the isotropy solutions 
provide also a good fit for the dispersion profile only (see 
e.g. the dashed lines in Fig. 2), they not correctly match 
the observed k. This produces a significantly worse to- 
tal x7dof (NGC 3379: 45/26; NGC4374: 35/40; NGC 
4494: 27/44) with respect to the best-fit in the Table 1, 
although still close to x 2 /dof <~ 1 mainly because of the 
large error bars. 

Finally, the best fit M/L+ in Table 1 are very similar 
to the values found for DM models (reported between 
brackets) in all cases, generally consistent with a Kroupa 
(2001) IMF. 

Looking at the f(R) parameters, in Fig. 3 we show 
the marginalized confidence contours of the main two 
potential parameters for the three galaxies. As also re- 
ported in Table 1, the 5 parameter has a mean value 
5 = —0.81 ± 0.07 which is inconsistent with the one pre- 
viously found for spiral galaxies (e.g. Sanders 1984, also 
shown in Fig. 3). On the contrary, S seems nicely corre- 
lated with the other potential parameter, L, as expected 
from Eq. 5. In the same figure the correlation is sup- 
ported by the tentative fit into the 5 — L plane (whether 
or not the spiral galaxy sample is included in the fit), 
although the sample is too small to drive any firm con- 
clusion. 

Interestingly, there seems to be a possible increasing 
trend of S with the orbital anisotropy: this is also shown 
in Fig. 3 where we have added the fiducial value obtained 
for the spiral sample (having assumed a reference /3 = — 1 
for fiducial tangential anisotropy for late-type systems, 
see e.g. Battaglia et al. 2005). This evidence leaves room 
for an interpretation of S and the physics of the galaxy 
collapse (e.g. the spherical infall model, Gunn & Gott 
1972; Gunn 1977). 

In fact, as discussed in §2, 5 is linked to 5i, which 
is an arbitrary function that comes out because the field 
equations in the post- Newtonian approximation, depend- 

10 For NGC 4374 only to be nicely fitted at all radii we needed to 
include some radial anisotropy in the very central regions, following 
the 0(r) distribution adopted in N+ll (see Eq. 5, whith best fit 
r a = 22.5). 
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Fig. 3. — Top: 1- and 2-<r confidence levels in the 8 — L space 
marginalized over M/L* and ft (see also Table 1). Spiral galaxy 
results from Sanders (1984) are shown as empty triangle with 
error bars. Solid (dashed) curve shows the tentative best-fit to 
the data including (excluding) the spiral galaxies and assuming a 
L <x tJ5/(1 + 5) correlation as expected from Eq. 5. Bottom: the 
anisotropy and the S parameters turn out to be correlated for the 
elliptical sample (full squares). This correlation seems to include 
also the spiral sample cumulatively shown as the empty triangle 
(here we have assumed = — 1.0 ± 0.5 as a fiducial value for spiral 
galaxies to draw a semi-quantitative trend across galaxy types). 

ing only on the radial coordinate. From a physical point 
of view, such a function could be related to second or- 
der effects related to anisotropics and non-homogeneities 
which could trigger the formation and the evolution of 
stellar systems. To take into account such a situa- 
tion, one should perform the post-Newtonian limit of 
the theory not only in the simple hypothesis of homoge- 
neous spherical symmetry (Schwarzschild solution) but 
also considering more realistic situation as Lemaitre- 
Tolman-Bondi solutions (see e.g. Herrera et al. 2010). 

5. DISCUSSION AND CONCLUSIONS 

There is a growing attention to alternative model to 
the ACDM paradigm as the latter is still suffering some 
discrepancies at the galaxy scales and, most importantly, 
is based on the assumption of the existence of two ingre- 
dients (DM and DE) whose nature is still unknown. 

Different attempts have tried to circumvent the prob- 
lem by introducing a modified dynamics, e.g. with the 
MOND theory (see Sanders & McGaugh 2002, Swaters et 
al. 2010, Cardone et al. 2010), but this seems still need- 
ing some DM at least to cluster scales which might be 
still consistent with the primordial nucleosynthesis (e.g. 
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TABLE 1 

Model parameters for the f(R) potential. 



Galaxy 


Mag (band) 


Res 


Ml Li, 


L 


5 


P 


X 2 /dof 


NGC3379 


-19.8(B) 


2.2 


6(7) 


6 


-0.75 


0.5(<0.8) 


14/25 


NGC4374 


-21.3(V) 


3.4 


6(6) 


24 


-0.88 


0.01(0.01) 


14/39 


NGC4494 


-20.5(B) 


6.1 


3(4) 


20 


-0.79 


0.5(0.5) 


18/43 



NOTES - Galaxy ID, total magnitude, effective radius and model parameters for 
the unified solution. DM-based estimates for M/L* and /3 (NGC 3379: DL+09; 
NGC 4374: N+ll; NGC 4494: N+09) are shown in parentheses for comparison. 
Ml are in solar units, R e g and L in kpc. Typical errors on M/L* are of the 
order of 0.2M/Lq and on /3 of 0.2 (see also Fig. 3). The small \ 2 values are 
mainly due to the large data error bars. 



via high energy neutrinos, Angus ct al. 2010) and does 
not provide an explanation for the DE. 

Lately /(-R)-gravity models have made their step out 
as a natural explanation for the two dark ingredients of 
the Universe assuming that they are related to the fact 
that gravitational interaction could present further de- 
grees of freedom whose dynamical effects emerge at large 
scales (Capozziello & Faraoni 2011). In this paper we 
have checked the Yukawa-like modification to the New- 
tonian potential obtained as post-Newtonian approxima- 
tion of /(i?)— gravity for which no particular choice of the 
Lagrangian has been provided, with the only assumption 
that f(R) is analytic function. 

We have used a combination of long-slit spectroscopy 
and planetary nebulae kinematics out to ~7 i? e ff in three 
systems (NGC 3379, NGC 4374, NGC 4494) for which 
ACDM models turned out to be fairly consistent with 
WMAP5 measurements (see N+ll for a discussion). 

Due to the small galaxy sample, the spirit of this anal- 
ysis has been to check whether 1) the modified potential 
introduced by the /(-R)-gravity allowed a fit to the galaxy 
kinematics comparable to the DM models; 2) the three 
galaxies returned a parameter 5 which is comparable with 
spiral galaxies (Sanders 1984). 

We have found that the modified potentials allow to 
nicely model the three galaxies with a distribution of 
the S parameters which turned out to be inconsistent 
with the results found in spiral systems. We have shown 
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some hints that S might be correlated with the galaxy 
anisotropy, /3, and the scale parameter, L, with elliptical 
and spiral galaxies following the same pattern. 

This evidence can have interesting implication on the 
ability of the theory to make predictions on the internal 
structure of the gravitating systems after their spherical 
collapse (e.g. Gunn 1977) which has to be confirmed on 
a larger galaxy sample which we expect to do in a near 
future. 

Despite of some simplifications on the model adopted 
(e.g. constant M/L and anisotropy across the galaxy) 
and the degeneracies between the model parameters, the 
results are very encouraging. The fit to the data is very 
good in all cases and both the stellar M/L (with Kroupa 
IMF generally favored) and orbital anisotropy turn out 
to be similar to the one estimated if a dark halo is con- 
sidered. 

Getting a modified gravity to work self-consistently for 
all gravitating systems in general, and all galaxy fami- 
lies in particular, is a very non-trivial challenge that has 
foiled other theories (e.g. MOND). 
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